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We report two analytical bounds for quantum error-correcting codes that do not have preexisting classical 
counterparts. Firstly the quantum Hamming and Singleton bounds are combined into a single tighter bound, and 
then the combined bound is further strengthened via the well-known Lloyd's theorem in classical coding theory, 
which claims that perfect codes, codes attaining the Hamming bound, do not exist if the Lloyd's polynomial has 
some non-integer zeros. Our bound characterizes quantitatively the improvement over the Hamming bound via 
the non-integerness of the zeros of the Lloyd's polynomial. In the case of 1 -error correcting codes our bound 
holds true for impure codes as well, which we conjecture to be always true, and for stabilizer codes there is a 
1-logical-qudit improvement for an infinite family of lengths. 



Introduction — During all kinds quantum informational 
processes there are uncontrollable noises due to inevitable in- 
teractions with ubiquitous environments. The theory of quan- 
tum error correction [lJ-|4J] provides us a powerful tool to pro- 
tect our precious quantum data from various noises by en- 
coding them into some special subspaces, called the quantum 
error-correcting codes (QECCs), which correct certain type 
of errors. Stabilizer formalism |5|-|9j] is one major method to 
find such kinds of subspaces, called stabilizer (additive) codes. 
There exist also some nonadditive codes IllOi 1 1111 with bet- 
ter p arameters for whose constructions a graphical approach 
1 12l Il3|] as well as an equivalent codeword stabilizer codes 
approach Il4ll has been developed recently. 

One of the fundamental tasks in the theory of QECCs is to 
find the optimal codes with as large as possible coding sub- 
spaces while correcting as many as possible errors and con- 
suming as less as possible resources. Obviously the optimality 
of the quantum codes, additive or nonadditive, is subjected to 
the tradeoffs among those parameters imposed by the princi- 
ples of quantum mechanics. The strongest bound so far is the 
quantum linear programming (qLP) bound II15[ 1 1611 which is 
unfortunately not analytical. So far all the analytic bounds for 



quantum codes have preexisting classical counterparts 1 17], 
for examples the quantum versions of Johnson bound 11811 . 
Griesmer bound 11911 . and most importantly, the Singleton 
bound (qSB) 0] and the Hamming bound (qHB) H. 

The qHB the qSB are two independent bounds and com- 
paratively the qHB is stronger for long codes and weaker for 
short codes than the qSB. For single error correcting codes the 
qHB is quite tight and can be attained [5, 6]. For instance in 
the case of qubits, i.e., two-level systems, all the optimal sta- 
bilizer codes saturating the qHB are constructed except a few 
families of lengths J2CJ Bill . Recently a family of nonaddi- 
tive codes attaining asymptotically the qHB has been reported 
1 2211 . For some families of lengths however the qHB admits 
strengthening II2(X|2U 12311 . which helps to identify some non- 
additive codes that outperform the optimal stabilizer codes. 

Here we will establish two analytical bounds without pre- 
existing classical counterparts for pure quantum codes. They 
are two strengthenings of the qHB with one arising from an 
interpolation with the qSB and the other one from a necessary 



condition for perfect codes attaining the Hamming bound. In 
the case of 1 -error-correcting codes we prove that the strength- 
ened qHB holds true for impure codes as well, which is con- 
jectured to be true for all distances. Though demonstrated in 
the case of pure quantum codes, both bounds apply to corre- 
sponding classical codes straightforwardly. 

Quantum Hamming-Singleton bound — In what follows our 
physical systems are qudits, p-level systems with p being fi- 
nite. In the Hilbert space of n qudits an A"-dimensional sub- 
space, whose projection is denoted by P, corrects a set of er- 
rors {Eu,} if and only if |4] (in an equivalent form) 

PE u tf a P = ^-Ti(P£ u £l)P. (1) 

If a set of errors is correctable then their linear combinations 
are also correctable so that we needs only to consider an op- 
erator basis, e.g., qudit Pauli errors {£^1 H. By a i-error we 
mean a Pauli error £ u acting nontrivially on < |w| = t qudits. 
As usual ((n, K, d)) p denotes a quantum code of length n (the 
number of the physical systems) of size K (the dimension of 
the coding subspace) and of distance d meaning that t-errors 
can be corrected, where t = [(d — l)/2j. 

The quantum code is said to be pure or non-degenerate if 
Tr(P£ u £ J) = whenever £^ ^ £ a , meaning that different 
errors send the original coding subspace into different orthog- 
onal subspaces. The existence of impure codes, where there 
are errors that we do not need to worry about, e.g., those errors 
that stabilize the coding subspace, makes an essential differ- 
ence between the classical and quantum codes. 

The quantum Hamming bound states that if a pure quantum 
code ((n, K, d)) p exists then (d = 2£ + 1 + a and a = 0, 1) 
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For odd distances (er = 0) the qHB can be established by a 
standard counting argument: there are iJ™ different i-errors 
and each on of such errors should take the coding subspace to 
a different orthogonal subspace. Though the above counting 
argument fails for impure codes, there are many positive ev- 
idences JUKHEltii] that the qHB holds true for impure 
codes as well. 



For even distances d = 2(t + 1) the qHB is usually taken 
as the same as that for d = 2i + 1. However a stronger bound 
as given in Eq.(f2| in the case of a = 1 can be read off readily 
from the classical coding theory, for which there exists also 
a simple counting argument. Consider the set of (t, l)-errors 
that includes all the /j-errors and all the (t + 1) -errors that act 
nontrivially on a fixed qudit. It is easy to count that there are 
Hiy of (t, l)-errors and the product of any two (t, l)-errors 
is a (2i + l)-error. Thanks to the error correction condition 
for pure codes each one (i, l)-error should bring the coding 
subspace to a different orthogonal subspace so that the qHB 
for even distance follows immediately. 

In contrast the quantum Singleton bound holds true both 
for pure and impure codes and it states that if a quantum code 
((n, K, d)) p exists then 

K < p n-2{d-i)_ (3) 

Though several proofs are known, including the original 
Knill-Lafflamme argument [4j] and the quantum LP bound 
16j,|9[], we provide here a simple argument which leads in what 
follows to a strengthening to qHB. Let us divide n qubits 
into 3 subsystems A, B, and C containing d — 1, d — 1, 
and c = n — 2(d — 1) qudits respectively. By introducing 
three reduced projections Pa = Tr^cP Pbc = Tr^P, and 
Pb = TtacP with P being the projection of the coding sub- 
space, we compute 
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KTrP 2 c > -TrP 2 , 
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in which the first equality follows due to the expansion Pa oc 
^ aCj4 Tr(P£ Q )£t since {£ a } a cA is a basis for subsystem 
A, the second equality is due to the error-correction condi- 
tion Eq.fQ}, and the inequality is due to the fact Tr(Pgc — 
Pb/p c ) 2 > 0. After interchanging subsystems A and B the 
above inequality must also hold so that the qSB follows. 

The first strengthening of the qHB comes from an interpo- 
lation of the qSB and qHB for pure codes. Let < e < t be 
an arbitrary integer and we divide n physical qudits into two 
subsystems A and B containing a = 2e and b = n — 2e qudits 
respectively. Denote 
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lu(ZB, \u\ <£ — e+cr 



where the summation is taken over all (t — e)-errors if a = 
while over all the (t — e, l)-errors if a = 1. From the counting 
arguments for the qHB it follows that TrP<. 
one hand we have 
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by noting that £^ J £ UJ i£ a is a (d — l)-error and the error cor- 
rection condition Eq.(Q]i for pure codes applies. On the other 
hand Tr(Tr A P e ) 2 > p" 6 (TrP e ) 2 so that we have 



Theorem 1 If a pure quantum code ((n,K, d)) p exists with 
d = 2t + 1 + a and a = 0, 1 then for all integers < e < t 



kh? : :<p\ Hi 



4eTrn-2e 
P n t-e,a- 



(7) 



Some remarks are in order, i) When e = t or e = 
the bound above coincides with the qSB or the qHB respec- 
tively so that our new bound provides a kind of interpolation 
and is called here as the quantum Hamming-Singleton bound 
(qHSB). ii) For a given length n the integer e can be chosen 
optimally according to 



e = t + 1 - 



P 
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with e = if n > tq 2 + 1 + a, i.e., it coincides with the 
qHB. As a result qHSB improves the qHB for short codes 
and in fact when t(p 2 — 2) > n — d > p 2 — 2 the qHSB is 
strictly stronger than both the qSB and qSB. We notice that 
n > 4(d — 1) so that there is no strengthening for qubits, i.e., 
p = 2. iii) Recalling that the codes attaining the qSB and qHB 
are called as maximum distance separable (MDS) codes and 
perfect codes respectively, the qHSB claims that there is no 
MDS code |90 if n > p 2 + d — 2 and no (pure) perfect code 
if n < d + t(p 2 — 2). In the next section the qHSB as well as 
the qHB is further strengthened via a necessary condition for 
the perfect codes to exist. 

Strengthened quantum Hamming bound — In the classical 
coding theory a crucial property of the perfect codes, codes 
attaining the Hamming bound, is described by the well-known 
Lloyd theorem 112511 which states that perfect codes of length 
n and distance d = 2t + 1 + a with a = 0, 1 do not exist if 
the Lloyd polynomial L™ a (x) = K™~ a (x — 1) has some 
non-integer zeros, where 



3=0 



(9) 



is the Krawtchouk polynomial of degree t. In what follows we 
shall demonstrate a quantitative version of the Lloyd theorem: 
how the nonintegerness of the zeros of the Lloyd polynomial 
gives rise to an improvement over the Hamming bound. 

Take an arbitrary integer < e < t and let Xj be the j-th 
root ( j = 1 , 2 , . . . , t — e) of the Lloyd's polynomial L^Z^a ( x ) 
of degree t — e. According to Ref. 12611 these roots are real and 
distinct with < Xj < n with their integer parts \xj\ being 
different integers. We introduce the following polynomial 



(x) 



n 
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l x j. 



(10) 



of degree 2(t — e). It is obvious that A"^(fc) > for in- 
teger k since the roots the polynomial A"^(a;) are pairwise 
consecutive integers. On the other hand A"^(xj) < for all 
1 < j < t — e and the equality happens for all j only when all 
the zeros Xj are integers. 



Theorem 2 i) If a pure QECC ((n, K, d)) p exists with d = 
2t + 1 + cr where er = 0, 1 then S™'^K < p n for all integers 
< e < t where 
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with i^^ e being the qHSB as defined in Eq.Q, x/s being the 

n — 2e / 
t — e,(7 \ 



zeros of Lloyd's polynomial L^_ e 2 ^.(x), and 



(*) - E 



[A'. 



n-2e- 
1 



-^-l)] 5 



J(p2-1)-1( 



_1 (n—2e- 



(12) 



being strictly positive, ii) The above strengthened qHSB is 
valid for impure codes as well for d = 3,4. 

Proof We shall postpone the proof of the first part to the 
Appendix since it is essentially classical and prove here the 
second part which is quantum mechanical since it involves 
impure codes. In the case of d = 3, 4, i.e., t = 1 and a = 0, 1 
we have e = and for simplicity we denote S^ = S™' a and 

A™(x) = A™' a (x). The Lloyd polynomial is linear and its 
single zero is determined by p 2 z a = (p 2 — l)(n — a) + 1. 
For convenience we denote <5 CT = 2 CT — l z a\, S a = 1 — 8 a and 
obviously (5 CT , <5er > 0. The strengthened qHSB is given by 
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According to the qLP bound if a code ((n, K, d)) p exists 
then there is a probability distribution {KA s /p n }™ =0 such 
that (Kg(x)) = A s for < s < 2 + a, where we have 
denoted (g(x)) = Sc Y^=o 9( S )A S for arbitrary g(x). Since 
{K"(x)} s < t provides a basis for the polynomial of degrees 
< t we can expand the polynomial A™ (a;) = (n — a;) 0- A" (a;) 
of degree 2 + cr by {A'"(:r)} s <3 and, when averaged with 
respect to the above probability distribution, we obtain 



<A2(z)> 



aoa + aigAj + ai a Ai + a 3a A 3 



(14) 



in which the coefficients are given by 

a 0(T = n^ilzal +S a - p 2 6 a S a ), 
a la = 2(n- l)"5 a +<7([z IJ \ + 8 a -p 2 S a 5 a ), 



<l2r, 



2(?i -2) a /p 2 + Aa8 a 



(Ho 



If we are able to show that ao CT A"(i) > aia- A™(0)(= n a ai a ) 
for all cr = 0, 1 and < i < 2+ct then the bound S[\ a K < p n 
follows immediately from the inequalities 

(K(x)) > ^-Y jL \l(j,)A i >^-Ya iu A i 
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Suppose at first n >p 2 + 2+cr. In this case we have [^o-J > 
p 2 and, since 4cviv < 1 and Ss,8 a > 0, [z a \ > p 2 S a S a + 



?>p 2 8 a /A so that ao a > 3n' T and ao CT > 4n cr 5 cr . Furthermore it 
is easy to see that A™(i) with fixed 1 < i < 2+cr as a function 
of L^o-J is increasing in the range [z a \ > 4 + a. By excluding 
a single case where p = 2, n = 7, and <i = 4 for which the 
bound is clear we can assume \z a \ > 4 + a so that A™(1) > 
(3 + a)/(5 + cr), A™(2) > (3 + cr)/10, and A™(3) > 1/5. If 
cr = it is now easy to check oooAq (i) > a^ for i = 1,2. 
In the case of cr = 1 it follows from 4n<5i(A"(l) — 1) > 
2(n — l)Si = ann — aoi (since n > 7) that aoiA™(l) > 
nan. In the mean time from aoi > rt(6/p 2 + 2c>i), 6A™(2) > 
2(n - 2), and 2(5^(2) > 4^ as long as n > 7 it follows 
that aoi A" (2) > 71021- Finally from aoi > in and n — 3 > 4 
it follows that o iA"(3) > 12n/5 > 6n/p 2 = nasi. 

Suppose now 4 + 2cr < n < q 2 + 1 + cr. In this case we 
have \_z a \ = p 2 S a = n — a — 1. As a result the coefficients 
can be evaluated as p 2 a,Q a = n a (n — a) [z a \, p 2 a\ a = (2 + 
a){n - l)(n - 2) CT , and p 2 a 2<7 = (2 + 4cr)(n - 2) CT . It is 
easy to check aoo-A"(i) > n° ' ai a for i = 1, 2, 3. In this case 
since p 2 a 0(T = n CT [^o-J ( W + 1) so that S^ CT = p 2 ( 2 + CT ) 
meaning that our strengthened bound coincides with the qSB 
K < "- 2 ( d - 1 ) in the case of d - 1 = 2 + a. Q.E.D. 

Some remarks are now in order, i) To obtain the corre- 
sponding strengthened Hamming bound for classical codes 
we have only to replace K by K/p n and then identify p 2 
with the number of the alphabet, ii) If the Lloyd's polynomial 
has at least one noninteger zero then |A™'^(a;j)| > so that 
S™ a :— S^' a > H™ a which yields the Lloyd theorem imme- 
diately, iii) When t = 2 andn = p 4 m((p 2 — l)m+ l) + 2 + cr 
with rn > all the zeros of the Lloyd's polynomial are inte- 
gers and there is no improvement. In the case of t > 3 (and 



since p > 4) 12711 there always exist non-integer zeros for 



Lloyd's polynomial so that we have always S" a > H™ a . iv) 
Numerical evidences show that the optimal value of the inte- 
ger e can be chosen as follows. Let {ajj}*_j be the zeros (in 
an increasing order) of L" cr (x) and let e = t — j with j being 
the greatest integer such that \_Xj\ < n — d + 2j. 

Let us now consider specially the stabilizer code whose size 
is K = p k for some integer k, which is the number of the 
logical qudits. The qHB and the strengthened qHB (e = 0) for 
pure stabilizer codes then read n — k > s p (n,d) > h p (n, d) 
where 



6<j/p 2 . (15) h p {n,d):=\log p H^], s p (n,d) := \log p S?J. (17) 



As long as S™ a > p h p( n > d '> we have s p (n, d) > h p (n, d) + 1 
since n, k, h p (n, d) are integers. In this case we say that the 
strengthened qHB has a 1-logical-qudit (1-l.q.) improvement 
over the qHB. Since p 2 S"™ = S^ 1 and p 2 H? = H 1 ^ 1 
a 1-l.q. improvement for the codes of length n and distance 
2£ + 1 is equivalent to a 1-l.q. improvement for the codes of 
length n+1 and distance 2t + 2. In general for a given d and n 
satisfying h p (n + 1, d) = h p (n, d) + 1 we may expect a 1-l.q. 
improvement unless all the zeros of the Lloyd polynomial are 
integers. 

Corollary s p (N^ a ,3 + a) = 2(?n+l + cr) while the qHB 
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TABLE I: Strengthened qHB S2(n,d) for binary pure stabilizer 
codes of length n < 128 that has a 1-logical-qudit improvement 
over the qHBs and meanwhile coincides with the qLP bound. 



for the same stabilizer code is 2 (to + a) + 1 where 
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with r being an integer satisfying 



< r < 
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(yi-4p 3 + 4p 4 - p 2 - (p - 1) 2 ) 



(18) 
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In the case of d = 4 and p = 2 the strengthened qHBs 
for lengths N^ x (to = 2,3,...) coincide with the qLP 
bounds as far as linear programming can be carried out. No- 
tably there is an infinite family binary stabilizer codes of 
lengths n a = (4 a — l)/3 with a > 3 for which we have 
^2(^0, 5) > h,2{n a , 5). And as far as numerical calculation is 
possible, these bounds also coincide with the qLP bounds. For 
larger distances the numerical results show that the strength- 
ened qHB for pure stabilizer codes also has a 1-1. q. improve- 
ment over the qHB as tabulated in Table I. There we have 
recorded all the strengthened qHBs for the codes of lengths 
< 128 and odd distances that have a 1-1. q. improvement over 
the qHBs and in the same time coincide with the qLP bounds. 
It seems that for relative long codes the strengthened qHBs 
with a 1-1. q. improvement will coincide with the qLP bounds. 

Conclusion and discussion — We have established two 
analytic bounds for the quantum error-correcting codes that 
do not have preexisting classical counterparts. Both bounds 
strengthen the quantum Hamming bound for pure codes and in 
the case of 1 -error-correcting codes they are valid also for im- 
pure codes as well. For stabilizer codes of a family of lengths 
the strengthened qHB has a one-logical-qudit improvement 
over the qHB. Nonadditive codes that outperform the stabi- 
lizer codes for these lengths may be expected and it is also 
interesting to find codes that attains the qHSB. Numerical 
evidences show that for relative large lengths the strength- 
ened qHB for stabilizer codes with a one-logical-qudit im- 
provement coincides with the qLP bound no matter how large 
the distance is. Furthermore the asymptotic behavior of the 



strengthened qHB deserves investigating. Finally we con- 
clude with a conjecture: 

Conjecture The strengthened qHSB KS™'* < p n with 
integer < e < t holds true for impure codes as well for 
arbitrary 
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Appendix: proof of Theorem 2 part i) — According to the 
quantum LP bound if a pure code ((n,K,d)) p exists then 
there exists a probability distribution {KA s /p n }™ =0 such that 
A = 1 and (K'"(x)) = for 1 < s < d where we have de- 
noted (g(x)) = |£ Z)"=o 9( s ) A s for an arbitrary g(x). 

Any polynomial g(x) of degree d — 1 has an expansion 
sl 3 -') = Si=o 9i^-?i x ) w i tn tne fi rst coefficient given by 
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n 



(20) 



where the inequality is due to A™'^(x) > for integer x. In 
order to evaluate the p— average on the left hand side of above 
equation we shall recall some properties of the Krawtchouk 
polynomials. 

First of all we have the Christoffel-Darboux formula 02611 

KHlO*£i(a0 " K™{x)K™ x {y) = 



m\x-ylA K™{x)K™(y) 



(t? - ^) t y V s { - ' 
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Due to the qLP bound we have (K™(x)) — for 1 < s < d 
from which it follows that (g(x)) = (g{x)) p as long as g(x) 
is a polynomial of degree d — 1. 

Take an arbitrary integer < e < t and denote r = 2e + cr 
for simplicity. We observe that (™~ :c ) A"'^(x) is a polynomial 
of degree d — 1 = 2t + a and thus we have 



Let Xj be the zeros of the Lloyd' polynomial L™_ 2e a (x) and 
xi ^ Xj if i ^ j. As a direct result of Eq.(f22b we have 



^ L^MWsZi:^) 



t (p 2 -i: 



,-i(»^) 



(23) 



((— )A- e (.)) p = (rr)A-(.))>^e) 



(21) Consequently the following polynomial of degree 2(t — e) 



,, , A „, e , , £?-£(*) ^ ^a-^-) ^ uizt.m^i-Ta^) 

J{ x ) - ^t,a\ x ) TTn - r 2^ „ . T t ( r -\ 2-^ 



ff 



t— e,0 



£j ayl^fo) ^ (pa _ iy 



1 /n— 2e— ct— 1\ 



(24) 



has a; = and x = Xj with j = 1, 2, ...,< — e as zeros so that we can write f(x) = x(p(x)L"_^(x) with ip{x) being some 
polynomial of degree t — e — 1. As it turns out, by denoting p™ = (p 2 — 1) 4 (") for convenience, 



t— e a n , e ( 



<(V)A?/ CT e (z)>, 



# 



^(W^W) + E A !^ L S-? } (<rK-i>))< 






P a t-efi 



V 



2(r+l) 



(25) 



Because of A"'^(xj) < for all j = 1, 2, . . . ,t and inequality Eq. (l2"Tl) the strengthened qHSB Eq. dTTb follows immediately. To 
carry out the calculations in Eq.(f25b we have used: i) The following two identities 



p X 



K^Z[-\x-l) 



-l)(n-r) 



I K s-i( x ) _ K" r (x) 



pZr 



(26) 



r 



for s > 1 and 



<(V) £ £*(*)>, = E<(V)«r- r (*)> P = ¥^ 27 ) 

which follow from two recurrence relations 



coefficients of u s v r in (1 — u) x (l + v+ (p 2 — l)u) n x in two 
different ways, and an identity 



Kr\x-l)=Y,K"{x). 



(30) 



p 2 x K?-\x-l) _ K?(x) _ K? +1 (x) 



(p 2 -l)r 



I't 



l>t 



Pt+i 



(28) 



ii) Due to the recurrence relations Eq. (12811 and Eq.(|29ll the 
polynomial x( n ~ x )L"Zea( x ) °f degree d can be expanded 
J-, ,,-.,-. by {K"(x) s > t — e\. Since the polynomial ip(x) of 

P { r ) K s (*) = HiJ( r-i ) K s+i( x )> (29) degree < t - e - 1 can be expanded by {K?(x) \ s < 



i=0 



t — e — 1}, the orthogonal relations (K"(x)K"(x)) f 



with the second recurrence relation obtained by computing the lead to (( n r x )f(x)) p = 0. 



